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ABSTRACT 



The specific impetus for this work was a conceptual 
design of a submarine using the toroid as the pressure hull. 
The designers could not find a ready body of knowledge to 
obtain scantlings for thier pressure hull. 

This work began with a review of efforts to solve 
complete toroidal structures. Several works were found which 
addressed general shells and extended into partial toroids, 
but. the solution of a complete toroid was not found to be a 
common exercise. Some of the these works are briefly 
reviewed. An attempt was then made to solve for the 
displacements in a thin walled circular toroid using the 
energy method. Several problems were identified associated 
with the structure geometry which make the solution for the 
complete toroid difficult. In addition, the functional used 
for the energy method needs to be more complex than the 
simple trigonometric or power series functionals used in this 
work. These two areas, geometry and functionals, are fertile 
areas for further study 
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CHAPTER 1 



INTRODUCTION 

The purpose of this work was to determine the response 
of a thin toroidal shell structure. The initial goal seemed 
achievable until the literature was reviewed. A simple and 
understandable solution was not available. The disparity in 
works was apparent in such basic elements as coordinate 
systems. The solutions were generally directed towards 
displacements or forces. The solution paths varied from 
complex algebra to finite difference method. The solutions 
were as varried as the number of individuals attempting them. 
It. was impossible to expect any simple correlation among the 
few people who attempted to address this topic. 

The result of this work is the identification of two 
areas where more time and effort is required. 

A. The impact of geometric curvature on the solution. 

B. The impact of the shell geometry on the assumed 
functional used in the energy method of the solution. 

To be very specific, these two areas have precluded this 
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effort from achieving the initial purpose of this- year long 



effort. 

A word needs to be addressed towards assumptions that 
are made. The assumptions are fairly standard for a shell 
being analysed in the linear elastic range. The assumptions 
are: 

1. The shell is made of isotropic and homogenious 
material which obeys Hook's Law. 

2. The thickness of the shell is constant. 

3. The thickness of the shell is small compared to 
the radii of curvature. 

4. A straight line normal to the middle surface 
before deformation remains straight and normal to the 
middle surface after deformation and retains its 
original length. 

In order to establish a consistent presentation, 
geometry will be discussed here and, unless otherwise noted, 
all representations of geometry will follow this coordinate 
system and the definitions. Solutions by other individuals 
will be translated into this reference system. 
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DEFINITIONS 



a : Ratio of R to r: a = - 

r 

Since R > r to form a complete toroid., a will always 
be greater than 1. 

C Circumference of the parallel circle; C = 2 n r* 

r : Radius of curvature in the 6 direction of the 

unloaded shell envelope, or radius of the circle 
which is rotated about the Z axis Cat distance R) to 
form the toroid. 

r : Radius to any point on the unloaded toroid measured 

perpendicularly to the Z axis. 

r* = R + r sin0 = r (a + sin0>. 

R Radius of rotation about the Z axis of the circle of 

radius r used to form the toroid. 

t : Thickness of the plate 

C> Angle of rotation about the Z axis measured counter 

clockwise from the positive X axis in the X-Y plane. 

0 Angle of rotation measured clockwise from the local 

perpendicular to the X-Y plane in the positive Z 
directon at distance R from the origin in any 
direction 0. See Figure 1. It sounds like a 
lot, but it is not. 
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Figure 1 
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Layout of Variables 
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P : Radius of curvature in the <p direction. 

p = — — + r 
sin© 

Coordinate System: The standard right-handed XYZ coordinate 
system is the global reference system. Other 
coordinates Ceg. © and <£) build on this basic 
reference system. 

Geometric Curvature: Curvature of the shell due solely to the 
geometry or curvature of the unloaded toroid. 

Load Curvature: Curvature of the shell caused by 

displacements in response to the applied load. This 
could also be looked upon as the change in curvature 
from the reference (unloaded) curved surface. 

Axis of Symmetry: Axis about which the small circle is 
rotated to form the toroid - the Z axis. 

Plane of Symmetry: Plane which bisects the toroid - the X-Y 
plane 

Meridian: Intersection of any plane containing the Z axis 

with the toroid. On either side of the Z axis, the 
intersection results in a circle (for the unloaded 
toroid) with © as the variable. (See Figure 2) 
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Parallel Circle: Intersection of any plane parallel to the 
X-Y plane with the toroid. (See Figure 2) 
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Figure 2 



Parallel Circle and Meridian 
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CHAPTER 2 



HISTORICAL VIEW AND THE PROBLEM 

Historical View: 

Upon starting to research the analysis of circular 
toroidal shells, the field was found to be not very large. 
Books on shell theory and analysis would neglect or only 
briefly mention the complete toroid. This could imply that 
the toroid was: 

< 1 ) a simple extension of general shell theory or 

<2) much more complex than general shell theory. 

The latter seems to be the case. 

Senjanovic tl] gives a brief history of shell theory 
in his book. His Russian pride shows in some instances 
<"... shortcomings were eleminated by representatives of the 
Russian School...”), but the relative youthfulness of the 
field is obvious. Initial work in what is now called shell 
theory was by H. Aron <1874) and A. Love <1888) a mere 
century ago. Senjanovic does cite Flugge - another work 
referenced here - as one of the developers of the classical 
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(.or western) theory of shells. Yet even Flugge devotes only four 
pages 12] to the toroid. The following works are illustrative of 
the background on toroidal shells. 

V. Flugge: Flugge starts with a generalized shell and 
establishes force equlibrium. From this equlibrium 
he solves differential equations to obtain solutions 
for forces and displacements. His discussion of the 
toroid points out some of the problems inherent in 
the full toroidal solution. 

L. Sobel: Sobel [3] was an understudy of Flugge. 

His doctorial dissertation starts with the same 
generalized shell theory but is slanted towards the 
buckling analysis of the toroid. His solution was 
worked on a computer and is difficult to follow. 

I. Senjanovic: Senjanovic starts with a force 
bal lance on a differential element in a localy 
orthogonal coordinate system, then uses several 
methods to solve them. Most notably, he employed 
complex variables and some computer solutions circa 
1970. 

S. Timoshenko: Timoshenko [43 starts with a plate 
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and develops into a generalized shell. His attention 



to the toroid is of the same order as Flugge. This 
writer found Timoshenko's general development of 
theory easier to follow than others and hence the 
references to Timoshenko in this area may be more 
frequent than to other authors listed here. 

E. Tsui: Tsui [5] starts with the constitutive 
relationships and develops them into forces. However 
for the toridal shape his solutions are on globally 
orientated forces and displacements. Tsui intended 
his work more as a handbook than a textbook. His 
solutions are tables of values from computer analysis 
of toroidal shells. He does provide however, a 
consistent development, though shorter and in less 
detail than some others, for various shells of 
revolution. 
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The Problem: 



The specific problem to be addressed in this thesis is 
the response of a toroidal shell to hydrostatic presure. 

More specifics on the loading will be addressed in chapter 
three. This is to set the stage for the problem itself 
and lay out the difficult areas. 

The goal of any structural solution is to take into 
account the loading, geometry of the structure, and the 
material in the structure and determine the response of the 
entire structure. The response is then compared to a desired 
standard or a set of failure criteria and adequacy of the 
dqsign is determined. This approach is generally slanted 
towards displacements of the structure since only through 
displacements are the applied loads distributed and 
equi 1 ibrated . 

In any structure, the various elements are in 
communication with one another. The beam under a simple 
tensile axial load is the easiest example, but a plate or 
shell must also satisfy this connectivity. Unfortunately, 
the complexity of this connectivity or communication among 
elements increases rapidly. In the axially loaded (tensile) 
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beam, only the axial displacement Csay the X direction) is 
important. Load the same beam laterally and now the X and Z 
directions are required to describe the beams response. 

Going to a plate now requires the X, Y, and Z directions. 

The shell now takes the plate out of the convenient XY plane 
for a reference. 

With appropriate coordinate system selection, some very 
practical shells can be handled easily. The cylinder and the 
sphere are classical examples. The transition to a toroid 
would - by simple implication - be very easy. Start with a 
cylinder of radius r and length of 2nR and bend it around 
upon itself. Discounting a few wrinkles on the inner 
sections of the cylinder, we now have a toroid. The wrinkles 
iiowever carry a much deeper implication. 

The structure's ability to distribute the applied load 
changes drastically. The complete toroid has no boundary 
conditions comparable to the cylinder's boundary conditions 
at X=0 and X=L. These are replaced by connectivity Cor 
continuity) requirements. The analysis of toroidal segments 
is facilited by the ability to insert boundary conditions. 

But even this facility can be of limited use if too much of 
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the structure is included in the analysis. The simple 
thin walled analysis on a cylinder for a force balance 
becomes an algebraic exercise of no small feat for the 
toroid. (This will be shown in chapter five.) 

It should be noted that all of the individuals above 
stated or implied some difficulty in obtaining a solution for 
a complete toroid. The most blatant was Tsui £6] when he 
states that in problems involving the complete toroid, the 
crowns ( O = 0°and 180° - see FIG 1) can not be solved 
due to singularities. Flugge £73 gives a solution for the 
forces but states that they cannot be used "... in the 
vicinity of the top and bottom circles” (Tsui's crowns) 
"without additional bending. . . ". He then provides a 
displacement expression containing terms which blow up at 
0 = 0° and 180°. 

With all of this going against an easy solution, two 

courses of action were undertaken: 

A. The singularities at the crowns could be avoided 

by using some expression which avoids the offensive 

term ( — — ) (see Chapter 6). This was avoided by 
sinO 

referencing the differential element to: 
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r = R + r sin © = r < a + sin 0 >' 

The biblical addage still holds: If thine eye offends 

thee, pluck it out. 

B. The method of solution selected would be one to 
avoid further problems as pointed out by Flugge [81. 
The energy method would avoid messy mathematical 
stumbling blocks in the solution, since the form of 
the solution is assumed going into the energy method. 
This was almost true as shall be pointed out latter. 
Additionally, this writer has not seen this attempted 
by others. Therefore, when the line forms to say 
that this method will not work, it should be a very 
short line indeed. 
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CHAPTER 3 



CONSTITUTIVE RELATIONSHIPS 

An expression, or series of expressions, relating the 
material motion to properties such as stress and strain (and 
further to forces and moments) is required. This follows the 
definition of' constitutive equation of the material as given 
by Fung 19], "...mathematical expression of the mechanical 
property of a material. " 

First we must define our displacements (see Figure 3). 
Three displacements will be used: 

u : Displacement in the positive © direction 
v : Displacement in the positive <p direction 
w : Displacement in the direction of the positive 
(ie. outward pointing) normal at the point of 
interest. 

Rotations will be defined as follows (see Figure 3): 

c >0 : Rotation of a unit normal about the tangent to 
the meridian in the positive u direction; that is. 



rotation about the positive © axis. 



Figure 3 



Dispalcements and Rotations 





Rotations 
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cjp : Rotation of a unit normal about the tangent to 
the parallel circle in the positive v direction; 
that is. about the local positive <p axis. 
a t . : Rotation about the positive w direction. 

The loading will impact greatly on the final outcome. 

For this problem, the loading will be assumed to be a 
constant hydrostatic pressure. This has several 
implications: 

A. Hydrostatic loading always presents a loading 
normal to the surface. 

B. Because of the global nature of the hydrostatic 
loading, any benifit to be garnered from symmetric 
conditions would be applicable to this problem. 

The second aspect of the loading., axisymmetric loading, 
results in the number of unknowns being reduced dramatically. 
Both Timoshenko [103 and Flugge till state that the total 
deformation can be represented by only two displacements; one 
normal to the surface and one along the meridian in the 
positive 0 direction. Flugge also states [121 that that 
stresses are independant of <p and in fact any derivatives 
with respect to 4> are zero. Timoshenko 1131 agrees, but 



21 



the derivative going to zero is implied., not stated as 



bluntly as in Flugge. 

To take this one step further, a lot of things will go 
to zero because of axisymmetric loading. Sheer stress and 
strain are cross coordinate terms, ie. 



Because of the type of loading described, v and o r are not 
required for the description of motion. v is not required 
because of the independence of the motion on <p, and o r is not 
required because it depends on a partial of 4>. All of this 
is the result of axisymmetric loading. 



Midplane strains are as follows: 

Cq : Is the strain in the positive 0 direction. Cq 
is made up of two parts Csee Figure 4): 

1 . Cq due to u 




But an\ 




At two points on the meridian; P and Q 



u 



p 



u 




e = £! = — — cu + — de - u) 



1 r d0 <50 



_ 1 6u 
r <50 



2. £q due to w 
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Figure 4 



Midplane Strains 
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Length of a segment is equal to the, radius 

times the angular rotation. The original 

length is r d©. At point P the radius 

increases by w to (r + w) so the new length 

would be Cr + w) d©. Similarly at point Q the 

radius increases but is now 

(r + w + — d©). The new length would be 
66 

Cr + w + ^d©)d©. Averaging the two lengths 
60 

at P and Q gives an average A1 = wd© + j 

Dividing by the original 1 = r d© gives: 

e = ^ tw + - — d©l Neglecting second 

r 2 6© 

order terms yields: E = f 

Total Cq is therefore: 

£ = 1 [*“ + „] 

6 r 66 

Cq : Is the strain in the positive <p direction 
(see Figure 4). Since the body is 
symmetrically loaded, all 4> f s have the same 
displacements. Therefore can be calculated 
by calculating the change in the circumference 
of the parallel circle at <p = constant. Since 
the circumference is related to the radius. 
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)J< 

the radius (r ) will be used to determine 



strain 



C = 2 n r 



* 



At point F“, displacements w and u contribute 



to the Ar as follows: 



Ar = u cos© + w sin© 



AC = 2 n Ar 



„ _ AC _ 1 



[u cos© + w sin©] 




ROTATIONS 

Oq is determined from two points, P and Q, on a meridian 
seperated by displacement u in the positive © direction Csee 
Figure 5). The unit normal vector at P must rotate 
through d© to go from P to Q due to displacement u, ie. 
r d© = u 



If w was constant from P to Q, no additional rotations would 
be required since the surface displaces uniformly through d©. 



required to maintain its normal relationship. This rotation 
will be: 



. • . d© 




If x 0.- then a rotation of the local normal at 0 will be 



'2 



£Hd©/(r d©)= - 
66 



1 6 \! 
v66 
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Figure S 



Rotation 




0 
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This makes the total 

<P 



_ 1 6w-. 

V r 60 

Because of the symmetric loading and no subsequent 
variations in the displacements with <p, no v displacements or 
— exists. However a can exist Csee Figure 6). 

Due to symmetry in loading, |w^| is constant in <p. 
However, for to change direction, an additional rotation 
perpendicular to must exist. This is the rotation 

ci0= - cos0 d 4> 

Figure 6 gives a graphical presentation of Oq. 

Timoshenko t 1 4 ] gives a similar result but achieves this as 



a side result of obtaining curvature. 



LOAD CURVATURE; 



Load curvatures are as follows: 

*0: Change in curvature along the positive 0 axis 
due to the applied load. 

Hp. Change in curvature along the positive 4> axis 
due to the applied load. 
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Figure 6 
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It is significant, to this writer at least, to seperate 



these curvatures from the geometric curvature defined in 

Chapter 1. The geometric curvature gives one an indication 

of how difficult the problem will be 1151 and defines the 

initial or unloaded reference surface for further work. The 

load curvature is the sole result of the structure’s response 

to the load applied and directly related to the displacements 

use to equilibrate the load. 

Load curvature along the 0 direction is the change in 

So, 

rotation about the © axis ( — *0 divided by the original 

60 

length (rd©). 



Timoshenko [16] 
similarly: 



*0 “ 



_ 1 



rd© 69 



— c 0rf ) d© 



Y 



Recalling: = - (u - — ) 

0 r 69 

„ = 1 2 <£" - 
w r 2 &e *e 2 

gives a description of this and is 



~ °© d©/(r* d©) 

_ cos © ( 6w-, 

r 2 (a + sin©) ^ 

Ve can now describe all the pertinent engineering 
parameters in terms of displacements. Adding to this list 
the general expression for planer stress, we get the 
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following: 






1-v 
_ E 



Ce e + ue <f? 



1-v 



5 4 ‘IS ♦ 111 ♦ 



°<P ~ 77^2 



— tu cos© + w sin©]) 

r(a + sin©) 

tu cos© + w sin©] + — t— + w3) 



i-f* rCa + sin©) r <5© 

This almost completes our constitutive relationships. 



Stress, strain, rotation, and local curvature can be 
expressed in terms of the displacements u and w. The last 
step is to proceed to normal forces and moments. This is as 



follows: 



h 

N = f 2 or dz ; V e ©, <p 

t _ V 

2 
h 

M = / 2 o z dz 

V -h 

2 
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CHAPTER 4 



EQUILIBRIUM OF FORCES 

The next aspect of the solution to be presented will be 
that of establishing equlibrium conditions for the 
differential element. The differential element is shown in 
Figure 7. Geometry rears it's ugly head in that the 
element needs two reference points to describe it's surface. 
The first is the axis of rotation of the shell in the © 
direction. The second is the origin of the global coordinate 
system. The first axis can be referenced to the global 
coordinate system by R and 0. 

Solutions for partial toroidal shapes all take advantage 
of the known boundary conditions on the edges. When a 
complete toroid is analyzed, the boundary conditions double 
back upon themselves and become internal vice external 
forces. In these instances, symmetric loading and 
deformations are invoked. Tsui 117] gives tabulated 
computer- results for both solutions and the differences are 
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Figure 7 



Differential Elevent 
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obvious. 



In the following discussion, and the associated figures, 
the following conventions will be used (see Figure 8>: 

F is a generalized force ie. a force or moment in 
the v direction (i =6, <£, n>. The units for F are 
t Force/unit length!. 

F is the force, ie. 

= [Force/unit length!*! length! ■ [Force! 

0 

N 0' Normal forces in the <^> direction acting on 

the side of the element on which <^> is constant. 

Sheer force acting on the 4 > « constant face of 
the element. Note: Due to symmetry of loading, 

= o. 

Moment to induce rotation in the <^> 
direction. The moment is defined as positive with 
respect to the right hand rule as applied about the 
local positive <^> direction. 

Each force will be addressed graphically with three 
isometric diagrams and a geometric breakdown to the 
appropriate F^ . The F^ for each force <N^, N^, . . ect.) can 

then be summed for a force balance. 
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Force 



Figure 



N 



e 



N 



4> 



P ( load) 

M, 



*e 



M 



<p 



8 

9 

10 

11 

12 
13 
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Figure 8 



Forces Defined 
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Figure P 



"© 




Contribution to F v doe to N 0 : 

> 

, > <5N > > 

F,: -N~ sin/3 + <N~ + — ^d6) cosa sin/3 
<t> O' 0 SO 

f 

> > > ' 

F~: -N~ cos/3 + CN a + — 2d©> cosa sin/3 

e © 1 0 <5© 

, , <sn' 

F : N d sina - <N~ + — 2d©> sina 

n © 0 6© 
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Figure 10 






V 




cx “ 0; a 



d6 



So. 



e. 



(3 = 0; ft «=d<£-o e - — 

&p 

Remember, — ■ 0 

exp 



Contribution to doe to N^: 

> > > > 

F -Njl + N * coso cosft 
<P <P <P 




*- & 



e 



F 0 : sina cosft - cosft sin a 



F n : N 0 sin P " N 0 sin ^ 
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Figure 11 



«e 




Contribution to doe to Q^: 



> > > > 



<p' 


Qq cosa sin/3 


+ cosa 


sin/3 




> 




> > 


> 


© : 


Qq sina 


cos/3 - 


Qq cosa 


sin/3 




> 




> > 


> 


l I 

n* 


Qq cosa 


cos/3 + 


Qq cosa 


cos/3 
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Figure 12 



P 




Contribution to F doe to Z: 

\ 

F^: - Z cosa sinf3 

F^: - Z cosf) sina 

F n : - Z cos/3 cosa 
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Figure 13 



«© 




Contribution to F doe to 



6 M 



F . : cos/3 - (M~ + — & d©) cos ft 

4> e u so 

, snL 

F~ “<M a + — ^d©> cos/3 sina 



e 



'© 



<5© 



F n : " n e sin/3 + <M e + 



6 M 

se 



©, 



© 



cosa 
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Figure 14 



M 









Contribution to doe to M^: 

, > . * * 
Fii M . cos/3 cosa - CM. + — ^ d©> cos/3 cosa 

<r <P & SO 

<5M. 



F~: - M . cos/3 sina + CM* + — & d©) cos/3 sina 
e <p <P 6$ 

> > . » > 
F : - M . cosa sin/3 + CM. + — E d©> cosa sin/3 

n 4> <P 
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Relationships between the N^and N must now be 
established. 



N$ = r* dO 



* N o > m Ar** 

N a + — - = (N A + — — ) <r* + d <f> 



<50 0 <50 <50 



<5N, 



= N A r*d<£ + n'c^I— > d0 + c.^t>T*dtp + < — §><££_>d* 
0 6 <50 <50 <50 <50 

, <5N’ * & * 

= CN^ + — — > r dtp, neglecting terms with— — as 



<50 



<50 



being relatively small 



Similarly: 



N . = N , r dO 
9 <P 



Q a = Q « r d 0 



o 



6Q, 



e, _ 



<*Q 



O a * 



Q a + < — -i = CQ A + — -> r dtp 
U 60 w <50 



M e = M 0 r dtp 
6 M, 



a m. 



M A + — - = CM A + — -) r* dtp 
0 60 0 60 

M . = Ml r dO 

<P •P 
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CHAPTER 5 



ENERGY METHOD 

Any attempted solution must maintain some sort of 
equlibrium. In the energy method, the equlibrium comes from 
minimizing the energy from the external load and the internal 
energy of the shell. For the external energy, a surface 
integral of the hydrostatic load and the normal displacements 
gives this value. The internal energy is determined by 
integrating stress times strain throughout the shell volume. 
The external energy is subtracted from the internal energy 
and the resulting functional is generally refered to as n, 
the total potential energy of the structure. With V 
representing the internal energy and V representing the 
external energy, the functional looks like: 

II = V - V 

If the value of II is at a minimum value, the structure 
wiil be in equlibrium. This is the principle of total 
potential energy 118]. This work will not address buckling 
or post buckling behavior, only the linear elastic region. 
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For this situation, equlibrium is achieved when: 

SB _ 6V_ _ 6V _ Q 
6q v 6q t 

where q^ is any generalized parameter used to determine V and 
V. 

Unfortunately, the parameters used to determine V and W 
are u and w - which is what we are attempting to solve. To 
get around this difficulity, a solution will be assumed. 

Then an energy balance will be conducted. 

The assumption of a solution is not without risk. 

First, the solution must meet known natural boundary 
conditions. Then the accuracy of the solution is limited by 
the closeness of t.he assumed solution to the actual Cand 
unknown) solution. 

In the case of the toroid, there are no "boundary 
conditions" as was pointed out in chapter four. In place of 
the boundary conditions, we have conditions of compatability 
to limit the solution. As stated earlier ([191 and 1201), 
symmetry of deformation can be assumed for the hydrostatic 
loading. The plane of symmetry is the XY plane. 
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To determine the internal energy we proceed as follows: 



V = / o c dV 

V v j v j 



V «. j 



e, 0 



£ CN £ + M * ) dS 

2S ij vj ij ij 



Recall that, due to the type of loading, 
cross terms Cv j ) equal zero. 



v - i S s < N e £ e - % c * ♦ M e V ds 



Eq 


1 : N 0 £ 0 - 


K 


< e © + ^ V e © 


; K 


= t n 

1-K- 2 




= 


K 


<£ 0 + 2 V £ 0 £^ + 


c p 






Simi larly : 










Eq 


2: 0 = 


K 


<.£^ + 2 P £^ £ 0 + 


e p 




Eq 


3: M e *© = 


[) 


C *6 + V *© 


> D 


_ E h 3 
12<l-v 2 ) 




= 


K 


C rz ) ( *© + 2 v *© 


** + 


*P 


Eq 


*+ = 


K 


C^) Cx 2 + 2 v 


*© + 


*P 



K <rz> 



Recall from chapter four: 



Eq 5,6: £ 0 = i (— + w) ; 



Cj - 



&e 



Eq 7,8: * 0 “ *a = 



^ r(a + sin©) 

2 v . ^ _ cos© 

2 J ? 



Cu cos© + w sin©) 



<u - ^) 



r“ <5© 6©^ ^ r^Ca + sin©) 6© 

The unknowns u and w still persist. The assumption will 

be made that u and w are functions of ©. A series solution 

will be assumed as follows: 
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w (6) = 2 w sinCm©}, limits of summation will be 

m 

addressed later 

u<0> = 2 u — (sinCmO)) 
m 60 

= 2 u m cos<m0) 

m 

This turns out to be not a very good assumption for the 
the series representation of the displacements. There are 
some inherent limitations; for instance sin<m0> is always 
zero at 0 = 0. While not in conflict with the conditions of 
compatabi lity , this is an additional constraint on the 
solution displacements. If the actual displacements are not 
zero at 0 = 0, then there is a guaranteed error in the 
assumed solution. The question to be answered, as with all 
assumed solutions, is how good is the answer obtained 
compared to the answer that is required. More on this 
subject in Chapter 6. 

For ease of notation, let: 

6 - sin(m0) 

m 

> 

0 = m cos(m0) 

m 

0 = -m 2 sinCm0> 

m 
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Therefore: 



Eq 9: w<6) = Z w <p 

m m 

Eq 10: u<©) = Z u <j>' 

m m 

Boundary conditions for this problem are actually 

compatabi lity conditions. At 6 ■ 90° and 270° the u<©> 

displacement should be zero and — <w<©>> should be zero. 

6 © 

= ^_<z w 0 ) = Z w — <0 ) « Z w (p 

m m rn td m m 

Since w are merely coefficients now in the series solution 

m 

they are unaffected by the derivative. To look at this in 
more detail: 

— = w cos<.6) + 2 cos<2©> + 3 cos<3©> + . . . 

60 1 2 3 

+ w n cos<n©> 

V i 

For conveince, let n = 10 for the trial solution. At © ■ 90° 
and 270°. cos(m©> = 0 for m = 1, 3, 5, . . and cos(m©> ■ ± 1 
for m = 2, 4. 6. . . . For the series to be equal to zero 

set = 0 for m = 2, 4, 6, . . . Similar reasoning gives 

u m * 0 for m = 2, 4, 6, . . . CThis explains why the 
functional for u was selected as the derivative of the 

m 

functional for w . > As shall be seen later, the method is 

m 7 

smart enough to realize this limitation.) 

The elemental area CdS) must be investigated for the 
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CR +' r sin©) d^ 



toroid. dS = dl 0 dl^ ; dl 0 <= r d© ; dl^ - 

= r d© (R + r sin©) d4> 

As a check: S « <R + r sin©) d© d tf> 

o o 

= 4 rt 2 R r. This agrees with the CRC 
Standard Mathematical Tables [213 for the surface area of a 
toroid. therefore dS is correct. 

The ultimate objective is to get * 0. In this case 

q = u and w . So for each specific u - that is u. - one 

v m m m k 

gets an expressison such that: 

6U _ 6V_ _ 6V_ m Q 

* u k 6u k 6u k 

Likewise a series of m expressions: 

SU _ 6V_ _ SV_ _ Q 

<5w, <5 w, 6w, 

k k k 

For a nominal m = 10, this yields 20 equations. Fortunately, 
u and w are 20 unknowns. 

rr» m 

Now the internal energy can be formated. Equations 9 
and 10 are substituted into equations 5 through 8. These 
equations are then inserted into equations 1 through 4. This 
is combined with the expression for dS to give V, the 
internal energy. To facilitate the future derivatives which 
will be required, the subscripts m and n will be used to 
indicate the seperate terms in a product, ie. : 
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= c a = v: ^ u <P + 2 w 0 ) ^ (Z u 0 + Z w 0 ) 

V7C7 I' m m m m P r» n r» n 

All ol' the above yields: 



V = n r 2 K 


/ 2n [-(Zu <p" + Zw 0 >i<Zu 0 “ + Zw 0 ) 

O m m m m r n n nn 


+ 

+ 


^(Zu 0 Zw 0 ) , , 1 ■ . _»(cos© Zu 0 + sin© Zw 0 ) 

I m m m~rri P C a + SXUOj n^n r» n 

R-rsTne) (COSe 2U nAn + Sine Z 


X 


rla i Sin0j CCOSe ZU n*n + sine Z "rAn 5 


+ 

+ 

+ 


C r? ) -T CZu m^ + s* *"> -i**" *" + s«A> 

J. 2 _ 2 mm mm ^ z ri ri r» ri 

r r 

<h^> 2^ CZu ^ + Zw 0^) cose <Zu 0* + £" n 0'> 

r r Ca + sin6) 

c'il) cos0 (Zu 0 * + Zw 0’> 

^ r 2 (a + sin©) 


X 


C , cos& (Zu 0’ + Zw 0*)] (a + sin©) d© 

r 2 (a + sin©) n n n n 


Next . 


carrying out the multiplications yields: 
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V * n r 2 K / 2n fi— IzXu u 00 + uw 00 + wu 00 

O m n m m mnmn mnmr 



+ W W 0 0) 

m n m n 



2 v 



r (a + sinG) 



• • > M 

IZXu u 00 cosG + u w 00 sinG 

m n m m mn m n 



+ w u 00 cosG + w w 00 sinG) 

mn m n mn m n 



— u u 0 0 cos G 

r 2 (.a + sir.G ) 2 m n m m 

> > 

+ u w 0 0 cosG sinG +■ w u 0 0 cosG sinG 

m r »v, m rs r rs 



m n m r, 



n m n 



+ w w 00 sin G) 

m n m n 



§ • t • 



+ C — )ZZ<.U U 00 “UW 00 - W U 00 

1 --, _ 4- m n m m mn m n mn m r 

zr 



+ W W 0 0) 

m n m n 

+ C--— -) 2 r cosG u 0 0 - uw 00 

12r 4 <a + sinG) mn rnm mn mn 
•• > *• > 

-WU 00 + WW 00) 

mn mn mn m n 

h 2 roc 2 fi > » > 

+ < n ) COS — ( U U 0 0 - U W 0 0 

2 mn m m mnmn 



I2r < a + sinG)' 
> > 



-wu 00 + w w 00)1<a+ sinG) dG 

mn m n mn r m n- J 



sv 

This is all well and good but actually - — is what we 

0 q v 

really need. Again, for this problem and ; and 

k = 1, 2, . . . 10. Applying term by term differentiation is 

straightforward but tricky and tedious. A couple of 
examples: 



— C w u 0 0 ) = w 0 0 

m n m n m m n 
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But observe: 



6u, 



(u„ U <P d> ) = u m d> rf>, 

m n m ~ r. m ™m “k 



U 0, 0 

n k n 



= 2 u <p 0. 

m m k 

Since m and n can both serve as a counter index, the two 
terms may be combined. The key to this is that both 
functionals C 0 ) have to be of the same order of the 

m 

derivative. Observe: 

U m K K * “n *i 

“m K * K 

Doing all of the above and collection like terms gives the 



S_ 

6u, 



C u u d> d> ) = 

TTi n m n 



following two expressions: 



The partial of V with respect to u k . 



6V _ 
6u, 



P U rr, 


/ 2n C2 <P 0, ■ 

O ^rrrk 


4 


2 cos 2 e 




Ca + sine) 2 


+ 


C^r)C2 0 m 0 k 


4 


2 cos 2 e 




Ca + sine) 2 


+Iw 

m 


/ 27t C2 
O k m 


+ 


> 

0 0 , cose) • 

^ m k 


+ c 


rZ )<_2 ^k^m 




2 cos 2 e 



2 v cos© + $’ $"') 

k m m 7 k 



sine) 






2 ^k 



’ 2 \> cos6 



. <t> 4>. )) Ca + sine) de 

2 T m T k 



2 



Ca + sine) 
2 cose sine 



C *k*\n Slne 



Ca + sine) 






2 ~k 



2 ^ ♦ «'>;> 



Ca + sine) 



Ca + sine) 

2 + ^e] 



rrr k ‘ 
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The partial oi’ V with respect to *# k . 



= n K flu / 2n (2 tp ' 4> 
6w L m ° 



2 P 



m k 



Ca + sin©) 



C 0 m 0 k sin© + 0 m 0 k cos©) 



2 cos© sin© j 



Ca + sin©) 



2 ^rrrk 



- c|X t )C 2 Ip tp 2 ^ cos© ( 0 ^ tp’ <p"'> 

17 ^m^k Ca + Sine) 'rn^k 



2 cos © 



Ca + sin ©) 2 m k 



ip ip )) Ca + sin©) d© 



*i. /*" (2 0 <p - - 4 - p s!? 6 - 0 A 

'" ° 1 m la ♦ sin©) k ” 



2 sin © 



Ca + sin©) 



2 KK 






2 v cos© 
Ca + sin©) 



C ip. ip + ^ <£, ) 

k m ~' 1 " 



rrr k ‘ 



2 cos 2 © 



<£. tp ))Ca + sin©) d©l 

km J 



Ca + sin ©) 2 k rn 
The external energy CV) should be looked at next. The 



energy is the product oi' the load - hydrostatic pressure 
CP) in this case which is always normal to the surface - 
times the displacement along the line of action of the force 



- wC©) - in this case. 



/ P w C©) dS 



_ J. 2 TT f 2Tt 0 rCR + r sin©) d© dip 

O o m ~m T 



Substitute R = r a 



= 2 n P r 2 Z w / 2Tr <p Ca + sin©) d© 

m O m 
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And the partials are: 



= 2 n P r 2 fV* <a + sin6) d6 
<5w, ° k 

k 

= o 

6u k 

The series expressions can now be rearranged: 

an _ <5V <5V 



<5q v 6q v ^ 
6V _ <SV . 

*q, *q, v 



o 



*V w k 



This can be represented in matrix format as follows 



Eq 11 



Q— = CD] <u> + [ E] <w> = 0 - — 



6u, 



<5u. 



XV 

Where [[>] represents the coefficients of u in 



6 u. 



(see page 9) 



IE] represents the coefficients of w k in 



<SV 

6u, 



< see page 9) 



<U> T = u = u , u , u , . u 

m 1 ' 2 ’ 3 ’ 1 0 



<W> T = « m = *V W 2 , W 3 , 



1 o 



Similarly in matrix format: 
Eq 12 



= [F] <u> + [G3 <w> = <P> = — 



6w, 



6 w. 



XV 

Where [FI represents the coefficients of u in - — 

<5w. 



(see page 10). 



[G1 represents the coefficients of w in 



6V 

<5w. 



(see page 10) 
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<p> T _ <5V _ 6V + 

<5w, <5w 

k 1 

Solving equation 11 

<u> 



6V 

6k, 



6w. 



2 lO 

for <u> gives: 



- [ D1 



- 1 



[El <w> 



Use this expression in equation 12: 

< [ FI t D1 " 1 [El + [G1 > <w> * <P> 



Solve for <w>: 



<w> = < t FI [ D 1 ” 1 [El + t G3 > 1 <P> 

If the solution for w and u are correct, then these 

mm 

can be used to reconstruct *K0) and u<©> and from these, 
using the constitutive relationships in chapter four, 
stresses in the 6 and <p directions can be calculated. Roark 
[221 gives expected values which can be used for comparison. 

Some justification for Roark's results should be 
provided since this is how the solution will be checked. He 
provides no specific reference for these answers in his book, 
but working backwards one can see that a simple thin walled 
approximation analysis was performed, similar to that of a 
cylinder. Roark provides: 




(See Figure 15) 




a + |sin0 
a + sin© 



> 
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Figure IS 



Thin Hailed Analysis 
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Op is a simple force balance in which the pressure times the 
internal cross sectional area equals the stress times the 
material area. 



PA - a. A 

<p m 

P (n r 2 ) » a* (.2 n r> t 

9 

o, = — 

* 2t 

A similar approach is used on Oq. The pressure acts on the 
area A. which in this case is <2 r)C2 n R). Oq acts at © and 
© + 180°. The thickness is constant at t but the 



circumference at © and © + 180° are different, 



= 2 n r = 2 rr r<a + sin©> 



The metal area at 0 is therefore A<©> = C<0) t. The force 



balance is: 



P A = o(0 ± ) CC0 t ) t +o C0 2 ) c<© 2 > t 

with ©„ = 0 and ©„ = 0 + 180°. Substitute Roark's o(©> into 
1 2 



the above equation gives: 

Pr a + l sine i 
P 4 rr r R = — C 2 1 



t a + sin© 



) 2 n rCa + sin©, > t 



p a + |sin© 

+ — C - -) 2 n rCa + sin© >t 

t a + sin©„ 2 

2 

Cancelling some terms shows the equality: 



2a «= 2a Q.E.D. 
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This must have been how Roark arrived at these relationships. 

This is how the solution will be checked. 

THE SOLUTION 

The method outlined above requires the integration of 
several quantities which are not readily found in a table of 
integrals. The integration was done using Simpsons Rule 
[23]. It was found, using known integrals, that 50 steps 
with in the interval 0 to 2 n gave good accuracy <5 to 6 
significant digits accuracy). A program was written to 
accomplish the integration and the subsequent matrix 
manipulations to solve for <u> and <w>. 

To check the output from the program, a spreadsheet was 
set up to handle: 



The spread sheet also calculated Oq, and o^. As a 



compare with the input pressure. 

THE RESULTS 

The results are disheartening. In the sample output 
provided, only one of the stress values used as a comparison 



w<.0> = I w 0 

m r i 



uC©> = i. u 0 

rvi i 



m m 



m m 




check, 0 ^ was averaged and used to calculate a pressure to 
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came close to the expected value. o ^ should have been 
constant and have a value of 5000 psi. It was not constant 
but It did manage to equlibrate most of the applied pressure 
<99.9%). This is pretty good for an approximation. The Oq 
is another story. 

Oq is too low and does not follow the solution in Roark 




The following are the default 
GEOMETRY 

Radius of rotation 
Radius of the circle 
Thickness of the shell 
MATERIAL 
Poisson's ratio 
: 3E+0? 

SOLUTION 
Number of steps in the Simpsons 
integral 
used in 

the series approximation 
Pressure 



a + sin© 

values in the program: 

: 8 inches 
: 2 inches 
: .02 inches 

: .3 Modulous of elasticity 



: 50 Number of terms to be 
: 10 

:-100 psi 
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m 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 



The program output for w and u are: 

m m 



w 

m 

-1 . 213117E-03 
6 . 043412E— 09 
-1 . 034239E— 04 
3. 250964E— 09 
-3. 933767E-05 
3. 328626E-09 
-1 . 16143E-05 
3 . 831 875E-09 
-2. 034983E-06 
2. 81 034E-09 



u 

m 

-1 . 137496E-03 
1. 491326E-09 
-1.151281E-05 
2. 026468E-10 
-1.579169E-06 
9. 227051E— 11 
-2. 382749E-07 
5. 9864E-11 
-2. 535156E-08 
2. 817199E-11 



Notice that the even numbered terms, which are to be set 
to zero as indicated earlier, are several orderes of 
magnitude below most of the odd values. 

The reconstructed output from the program is as 
follows: 
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The spreadsheet out put to check the solution is: 





P = 


-100. 00 


psi 


-99.90 < 


Sin Press 


Cpsi) 




V = 


0. 30 












E 

1-v 2 


3. 3E+07 


E = 


3E+07 








k - 


8. 00 


in a = 


4.00 








r = 


2. 00 


in 










t = 


0. 02 


in 








© 


w 


6m 

60 


u 


<5u 

60 


°© 




<°> 


C in> 


( raT } 


( in) 


.in . 
^racr 


Cpsi) 


<psi> 


0 


0. 0E+00 


-1.8E-03 


-1 . 2E-03 


0.0E+00 


-1461 


-4870 


15 


-4 . 4E-04 


-1 . 4E-03 


-1 . IE-03 


4. 2E-04 


-1710 


-4738 


30 


-7. 2E-04 


-8. IE-04 


-9. 8E-04 


6.8E-04 


-1946 


-4607 


45 


-9. 0E-04 


-5. 7E-04 


-7. 8E-04 


8. 4E-04 


-2123 


-4405 


60 


-1 . 0E-03 


-4 . 2E-04 


-5. 4E-04 


9. 6E-04 


-2258 


-4248 


75 


-1. IE-03 


-2. 2E-04 


-2. 8E-04 


1 . 0E-03 


-2345 


-4159 


90 


-1 . IE-03 


-2. 3E-08 


-2. 8E-08 


1. IE-03 


-2375 


-4131 


105 


-1 . IE-03 


2. 2E-04 


2. 8E-04 


1.0E-03 


-2345 


-4159 


120 


-1 . 0E-03 


4 . 2E-04 


5. 4E-04 


9. 6E-04 


-2258 


-4248 


135 


-9. 0E-04 


5. 7E-04 


7. 8E-04 


8. 4E-04 


-2123 


-4405 


150 


-7. 2E-04 


8. IE-04 


9. 8E-04 


6. 8E-04 


-1946 


-4607 


165 


-4. 4E-04 


1 . 4E-03 


1. IE-03 


4. 2E-04 


-1710 


-4738 


180 


-9. 6E-08 


1 . 8E-03 


1 . 2E-03 


9. 2E-08 


-1461 


-4870 


195 


4 . 4 E— 04 


1 . 4 E— 03 


1 . 1 E— 03 


-4. 2E-04 


-1269 


-5190 


210 


7. 2E-04 


8. IE-04 


9. 8E-04 


-6. 8E-04 


-1085 


-5498 


225 


9. 0E-04 


5. 7E-04 


7. 8E-04 


-8. 4E-04 


-894 


-5654 


240 


1 . 0E-03 


4. 2E-04 


5. 4E-04 


-9. 6E-04 


-747 


-5769 


255 


1 . IE-03 


2. 2E-04 


2. 8E-04 


-1 . 0E-03 


-659 


-5854 


270 


1 . IE-03 


6. 8E-08 


8. 4E-08 


-1. IE-03 


-630 


-5886 


285 


1 . IE-03 


-2. 2E-04 


-2. 8E-04 


-1. 0E-03 


-659 


-5854 


300 


1 . 0E-03 


-4. 2E-04 


-5. 4E-04 


-9. 6E-04 


-747 


-5769 


315 


9. 0E-04 


-5. 7E-04 


-7. 8E-04 


-8. 4E-04 


-894 


-5654 


330 


7. 2E-04 


-8. IE-04 


-9. 8E-04 


-6. 8E-04 


-1085 


-5498 


345 


4. 4E-04 


-1 . 4E-03 


-1. IE-03 


-4. 2E-04 


-1269 


-5190 


360 


1 . 9E-07 


-1 . 8E-03 


-1 . 2E-03 


-1 . 8E-07 


-1461 


-4870 



Note: Due to the precision of the machine being used, the 
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forced compatability 
forced conditions- at 
than the rest of the 



conditions appear to be non-zero. The 
e = 90° and 270° are significantly less 
solution and therefore satisfactory. 
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CHAPTER 6 



CONCLUSIONS 

The apparent simplicity of the toroidal shell belies the 
complexity of the mathematics to analyse this shell 
structure. The analysis of simpler shells Cplates, 
cylinders, and spheres) has the advantage of the geometric 
curvature remaining either constant or at lease maintaining 
its sign positive or negative. The geometric curvature of 
the toroid changes from positive to negative as 0 goes from 0 
to 2 it. 

The impact in the change of sign in the radius 
of curvature about the 2 axis has been seen in the results 
of this work and in applying other's solutions. This impact 
is emperical and could be the subject for further 
investigation. The first blatant statement was found in Tsui 
[24] where he states, "...the relevent differential equations 
to assume singular solutions, and consequently, problems 
involving the crowns such as a complete toroid cannot be 
solved. " The culpret was writing the differential equation 
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in terms of (Tsui's nomenclature): 



r = b (1 + — — ) 
sin© 

<This papers nomenclature): 

p = r (1 + (^/sin©) 

To avoid the — term, Tsui shifts variables - both 

sin© 

dependant and independant - and provides influence 
coefficient matrices. (These are outputs from numerical 
analysis. ) 

As mentioned earlier, Tsui addresses the significant 
differences between the partial toroidal shape and the 
complete toroid. Others may have addressed partial toroidal 
shapes tangentially but Tsui once again is very blunt in his 
statements. Again remember Tsui is intending his work as a 
practical handbook. His intended audiance wants practical 
answers to real problems. Others, Flugge [251 
specifically, attempt to give solutions but he adds, "...this 
solution cannot be realized. .. because it again leads to an 
incompatability of deformations. . . **. 

CURVATURE 

In attempting to point the finger at why the solution 
is so difficult, the answer comes back to curvature. This 
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work has seperated curvature into geometric (.unloaded) and 

loaded curvature. This writer feels that the geometric 

curvature is the culpret. Notice the expression for the 

radius of curvature in the <t> direction given in chapter one: 

R 



sin© 



+ r 



For the range 0° < © < 180°; sin© > 0 and p > 0. Conversely, 
for 180° < © < 360°; sin© < 0 and p < 0. At © ■ 0° and 180° 
sin© = 0 and p = ± ». Two things were observed. 

A. In previous solutions obtained in this study- ie. 
worse than the one presented - the stress outputs for 
the two regions differed markedly between 0° < © < 



180° and 



180° < © < 360°. This could have 



been due to a lot of reasons, but this was also 
coincident with . . . 

B. Flugge’s solution for the u displacement [26] 
was being investigated. This tended to blow up from 
© = 180° to 360°/0° due to a term lnCtan^). Other 
problem terms existed <cot©> which agrivated the 
solution at the points © = 0° and 180°. In just 
"playing” with the solution, the points of curvature 



change kept comming back as critical points requiring 
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further attention. 



Something happened at 0 «= 0° and 180°. In a cylinder 
these two points are no trouble - but the radii of curvature 
are «> and r - and are constant! In a sphere these two points 
are no problem - but the radius of curvature is r in all 
directions - and constant! In the toroid the radii of 
curvature at both points are p = ± oo and r. The key 
difference is that p is changing from p > 0 to p < 0 and at 
6 = 0° and 180°. Cor just happens to be passing through 
zero. 

It was mentioned earlier in Chapter 4 that it required 
two reference points to define the toroid. Some individuals 
may take exception to that statement, but here is where it 
comes into play. Another way to look at p, the radius of 
curvature in the </» direction, is that p is the length of the 
arm joining the point in question Csay point P) and the axis 
of symmetry Cthe Z axis in our case) and coincident with the 
local unit normal vector at point P. The angle 0 can be 
measured either at the axis of rotation of the small circle 
forming the shell, or at the intersection of the line p and 
the axis of symmetry for the entire toroid. Now when 0 ■ 50° 
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lor instance, the geometry looks like that shown In Figure 
16. As © decreases toward 0, tRe end of p must travel 
farther and farther down the axis of symmetry. All this time 
p is getting larger and is getting smaller. When © gets 
to zero plus <0 + ) the situation is one of approximations. p 
is the hypotenuse of the right triangle with the right-angle 

p 

at the origin. Sin© = £ from simple trigonometry, but as p 

h' 

goes to oo, sin© goes to zero and hence © must go to zero 
also. 

Continuing © around to 0~, the closest point on the axis 
of symmetry in line with the normal vector is now at +oo. 

This now gives p a direction opposite to that of the local 
normal vector, hence a minus sign, and as © goes further 
negative, p is now smaller in magnitude. (Similar arguments 
can be made for the bottom crown where © * 180°) The sole 
reason for the ± «> trip in p, and hence the sign change in 
is the offset R. Somehow this must be related to the 
difficulty in solving the complete toroid. 

GEOMETRY 

The other subtility of the toroid is the way the 
mathematical representation of the surface impacted on the 
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attempt- to integrate for energy over the toroid surface. 

This work specifically avoided some problems encountered by 

other investigators, significantly the expression for the 

curvature in the <j> direction. The early development in this 

work specifically avoided the term C — — ) , selecting instead 

sin© 

an expression containing the term Ca + sin©). Force balance 
and constitutive equations were developed using the latter 
expression. When the constitutive equations were plugged 
into the energy expression CV), no singularities existed. 

This was good news. 

Unfortunately, the Ca + sin©) term also appeared in the 
expression lor the differential surface area CdS). The bad 
news was that, many of the elements in V contained 
combinations of trigonometric functions not readily found in 
the CRC Math Tables [271. This is what forced the numerical 
integration of the elements in V. 

This heavy dependence of the trigonometric functions in 
describing the problem has an astounding impact on the 
selection of an assumed functional for the energy method 
solution. Any combination of CA © sin©) or C© A sin©) goes to 
zero when integrating from 0 to 2n. Remember that the dS 
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term contains a sinO term so all integrals contain at least 
one sin© term. 

It should be pointed out that the solution to a 
symmetric problem can be obtained by solving only one of the 
unique parts of the problem. This was not done in this case 
because it was hoped that the program being developed would 
be more general in nature, actually leading to orthotropic 
analysis of a stiffened toroid. 

The only redeeming quality of the integrals also comes 
from geometry. Geometry terms containing cos©, cos 2 ©, sin©, 
sin 2 ©, and Ca + sin©> and thier products appear in numerators 
and denominators through out the integrals. Mister Simpson 
and the digital computer were allowed to seperate the wheat 
from the chaff. 

An attempt was made to develope a non-trigonometric 
functional which would satisfy the conditions of 

connectivity stated earlier and to meet connectivity at 
© = 0°/360°. The conditions to be met were: 
f<0> = f <360) 
f'cO) = f ’ <360) 
f C 90) = 0 
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f <90) - 0 



f <270) * 0 
f <270) = 0 

This resulted in eight conditions. A functional of the 
following form was generated: 

f<m0> - A + B<m0) + C<m0) 2 + D<m0> 3 

+ E<m0)* + F<m0> 5 + G<m0> 6 + H<m0> 7 
A was arbitrarily set to 1 <to force a non zero solution at 
0 = 0°). Unfortunately, by the time f<mq> was processed, 
the functional looked the same for all values of m. The 
coefficients all changed as m changed, but the shape of the 
fuhctional remained unchanged from one value of m to the 
next because of the eight conditions imposed upon the 
functional. This lost the varietj' of the solution that one 
expects to obtain using a varying m. See Figure 17. 

It is obvious that more complex functionals are needed. 
This work included an initial evaluation of several 
functions for this task. A considerable amount of time was 
put into investigating the hypergeometric functions, 
particularly the Legendre functions. In addition to serving 
as a functional for u and w, it was also hoped that this 
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function would satisfy the form of the differential equation 



for the toroid and hence yield an exact solution. <It was 
partially the anticipation of using such a complex function 
that drove the layout of the program (see Appendix) to have a 
seperate generating section for the functional.) The 
Legendre function, or any other hypergeometric function, was 
not used due to compatability conditions being harder to 
achieve. The next step in this progression would be the 
assumption of a combined trigonometric function for u and w, 
Cie. A sin<m0) + B cosCm©)). 

m m 
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Figure 1 6 
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Figure 17 



Generated Functional 
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CLOSING 



In Chapter 1 it was stated that the geometric curvature 

and the impact of shell geometry on the assumed functional 
were the two problem areas in this thesis. Of the two, the 
former is the more significant in this writers opinion. The 
analysis behind Figure 15 is crude, but it is an attempt to 
represent in graphics and mathematics what the shell is 
doing. Being unable to do this is what makes the solution so 
difficult. 

In a dynamics class, one student asked the professor if 
another mathematical technique could be used to analyse a 
spring. The professor thought briefly, then answered, "You 
can do all the math you want but you have to think like a 
spring. The spring knows what it is doing. " This writer 
feels a little short on thinking like a toroid, but also 
feels in good company. 
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APPENDIX 



THE PROGRAM 

This program is to set up the geometry, the energy 
expression, Simpson’s multipliers, and perform the matrix 
manipulations to solve for the displacements Cu> and Cw>. 



50 CLS: PRINT: PRINT: INPUT "WHAT OUTPUT FILE NAME ";AS 
60 CLS 

1000 NN= 50 ♦This is the number of steps to be used in the 

Simpson's integral 

1010 MM=1 0 ♦This is the number of terms in the series for 



1020 

1030 

1040 

1050 

1060 

1070 

1080 

1090 

1100 

1110 

1120 

1130 

1140 

1150 

1160 

1170 

1180 

1190 

1200 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 



displacements u ans w. 

NU=.3 ♦Poisson's ratio 

R=2 ♦R from the text. 

A=4 ♦a from the text. 

H=.02 *t from the text. 

E=3E+07 ♦Young’s modulus 
P=-100 ♦pressure in psi. 

CLS: PRINT: PRINT: PRINT 

PRINT “The following are the default values in the 
program: " 

PRINT: PRINT 
PRINT “GEOMETRY" 

PRINT "Radius of rotation 
PRINT "Radius of the circle 
PRINT "Thickness of the shell 
PRINT: PRINT 
PRINT "MATERIAL- 
PRINT "Poisson's ratio 
PRINT "Modulous of elasticity 
PRINT: PRINT 
PRINT "SOLUTION- 
PRINT "Number of steps in the Simpsons" 

PRINT " integral :";NN 

PRINT "Number of terms to be used in" 

PRINT " the series approximation :";MM 

PRINT "Pressure :”;P;" psi" 

PRINT: PRINT 

INPUT "Do you wish to change any of these CY/N): ":SS 

IF S$="N" THEN 1520 

CLS 

PRINT: PRINT 
PRINT "GEOMETRY" 

INPUT "Radius of rotation Cinches) :";RR 

INPUT "Radius of the circle Cinches) :";R 

A=RR/R 

IF A > 1 THEN 1380 

CLS: PRINT"RADIUS OF ROTATION MUST BE GREATER THAN THE 



” : A»R: " inches" 
” : R: " inches" 

"; H; " inches" 




RADIUS" 

1370 PRINT "OF THE CIRCLE. ": GOTO 1300 
1380 INPUT "Thickness of the shell Cinches) :";H 
1390 PRINT: PRINT 
1400 PRINT "MATERIAL- 

1410 INPUT "Poisson's ratio :";NU 

1420 INPUT "Modulus of elasticity :":E 

1430 PRINT: PRINT1440 PRINT "SOLUTION" 

1450 PRINT "Number of steps in the Simpsons" 
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1460 INPUT " integral 

1470 PRINT "Number of terms to be used in" 
1480 INPUT " the series approximation 
1490 INPUT "Pressure Cpsi) 

^ Rnn > 



1500 ’ 
1510 •’ 



integral 




1520 ’THIS SECTION IS TO SET UP THE PARAMETERS FOR THE 
INTEGRATION SECTION 

1530 ’ 

1540 DT=6.283185307#/NN 
1550 I=NN+1 

1560 DIM PC 1. 3, MM) ^Defines the functions <p. 

1570 DIM JCI) *Simpson’s multiplier 

Arrays MD, ME. MF, and MG correspond to the matrixes in 
the partial of V with respect to u k and « k . MP is the matrix 

representing the energy due to the pressure term. MS is a 
scratch matrix. M*~ is a series of matrixes used in the 
matrix manipulation subroutines. MV and MU are the answers. 

1580 DIM MDCMM . MM) : DIM MECMM, MM) : DIM MFCMM, MM) 

1590 DIM MGCMM. MM): DIM MSCMM, MM) 

1600 DIM MX-CMM . MM) : DIM MY#CMM, MM) : DIM MZ#CMM, MM) 



1610 DIM MPCMM) : DIM MVCMM) : DIM MUCMM) 

1620 CLS: PRINT: PRINT: PRINT "GENERATING FUNCTIONS" 
1630 FOR I = 1 TO NN+1 
1640 T=C 1-1 )*DT 
1650 FOR M = 1 TO MM 

A D/T A M \ M 7 ? Tk ; „ 4^. 4 _ 4^. . 



1690 NEXT M 
1700 JCI) =4 

1710 IF ( 1/2-INTC 1/2) ) >0 THEN JC I )=2 

1720 NEXT I 

1730 J(1)=1:J(NN+1)=1 

Integrate using Simpson’s Rule 

2000 ’INTEGRATE FROM 0 TO 2*PI 
2010 IC=3 . 14159*CE*H/Cl-NU"2))*DT/3 
2020 BC=H'2/C6*R'3) 

2030 PC=-2*3. 14159*P*R'2*DT/3 
2040 FOR I = 1 TO NN+1 
2050 T=CI-1)*DT 
2060 S=SIN(T) 

2070 C=COS(T) 

2080 CLS: PRINT: PRINT: PRINT"INTEGRATING ON STEP ":I 
2090 PRINT " ANGLE ": T*57. 296; " DEGREES" 

2100 AA=A+S21 02 HH=H' 2/C12*R~2) 

2104 CC=C/AA 
2106 CA=2*NU*CC 
2108 CB=CC' 2 
2110 CD=2*NU > J < S/AA 
2112 AX=2*S*C/AA' 2 



: DIM MI^CMM.MM) 



1660 PCI.1 ,M)=SINCM*T) 

1670 PC 1 . 2 , M)=M*C0SC M*T) 

168-0 PCI.3,M)=-M'2*SINCM*T) 



’This is the function 
’This is the first 
derivative 
’This is the second 
derivative 
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2130 

2140 

2150 

2160 

2170 

2180 

2190 



2200 



2210 



2220 



2230 

2240 

2250 

2260 

2270 

2280 

2290 

2300 

2310 

2320 

2330 

2340 

2350 

2360 



FOE K = 1 TO MM 
MP(R)=MP(R)-J(I)*P(I, 1 , R)*AA 
SQ=1/AA' 2 
TR=2*NU/AA 

BX=H ' 3*( 1 +2*NU*C/AA+( C/AA) ~ 2 ) /( 1 2*R~ 2 ) 

FOR M = 1 TO MM 

MDCK, M)=MD(R, M)+J( I)>*((1+HH)*(2*P(1, 3, K)*P(1, 3, M) 

+CA*Pa.3,K)*P<:i,2,M$ 

+CB*2*P( I 2,D#P(I,2,M)))#AA 
MECR, M)=ME(R, M)+J( I)*(2*f*( t 1 M)*f*( i, 3, K) 

+CA*(Pa,3 i E$*f > <:i.i .m$*§ 
+Pa,l,M>*£(I,2,0*fc) 
+AX*PCI.2.K)*P(i,l JD 
-HH*(2*P( 1, 3, M)*p(i , 3 K) 

+CA*(f > ( t 3 . 2 , M) 

+P( I . 3 . m5*^ci . 2 . k5) 
+2*CB*P(I,2.k5*P(I.2,M)))*AA 
MF(R. M)=MF(R. M)+J(I)>K2*P(I,3, M)*P(I. 1,K$ 

+2*NU*(P( 1 . 1 . K5 *PC 1 . 3 . M)*S 
+PCI,2.M)*PCI,1,K)«<C)/AA 
+AX*P(I.2,M)*P(I,1 . K) 

-HH*(2*P( t . 3 . M)*P( 1 . 3 . K) 

+CA*(P( I, 3,M)*P( I, 2 , R) 
+P(I.3 i K$*P(I,2 i M$) 

+2*CB*P( 1,2. M$*P(I, 2, K) ) >*AA 
MG(R.M)=MGCR,M)+J(I)*(2*PCI,1,M5*P(I i ,k5 

+2*CD*P(I 1 , M)*P( 1.1 ,K) 

+2*P(I.l .fl)>KPCI.l,K)*S"2/AA"2 
+HH*C2*P( 1,3 . M)*P(I . 3 , K) 

+CA*(f > ( 1, 3,M)*P(I . 2, R) 

+PCI . 3 . K.$*£( 1 . 2, M5) 

+2*CB*P( 1,2, M$*k I , 2 , K) ) )*AA 

NEXT M 
NEXT R 
NEXT I 



FOR K = 1 TO MM 
MP(R)=PC*MP(R) 

FOR M = 1 TO MM 
MD(R,M)=IC*MDCK,M) 
MECR, M)=IC*ME(R, M) 
MFCR.M)=IC*MF(R.M) 
MG(R.M)=IC*MG(R,M) 
NEXT M 
NEXT R 



3000 ’SOLVE THE MATRICES FOR w AND u 

m m 

3010 CLS: PRINT: PRINT: PRINT ••SOLVING" 

3020 ’ 

3030 FOR R = 1 TO MM: FOR M = 1 TO MM: MX#(K, M)=MDCR, M) : NEXT 
M: NEXT R 

3040 GOSUB 10000 ’ ***** Invert MD **** 

3050 FOR R = 1 TO MM: FOR M = 1 TO MM 
3060 MY#(R, M)=ME(R. M) : MX^CK, M)=-MI#(K, M) 

3070 NEXT M: NEXT R 

3080 GOSUB 11000 ’ **** (-MD'-1)*ME **** 

3090 FOR R = 1 TO MM: FOR M = 1 TO MM 
3100 MIXR. M)=MZ^(R. M) : MXXR. M)=MFCR, M) 

3110 MSCR. M)=MZ^(R. M) » **** Saving this for Urn 

solution **** 

3120 NEXT M: NEXT R 
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3130 GOSUB 11000 ’ **** MF*C-MD'-1 )*ME »*o>okh< 

3140 FOR K = 1 TO MM: FOR M = 1 TO MM 
3150 MY#CK,N)=MGCK,M): MX^CK, M)=MZ#CK, M) 

3160 NEXT M: NEXT K 

3170 GOSUB 12000 ' **** MG - MF*C-MD~-1 )*ME •Mow 

3180 FOR K = 1 TO MM: FOR M = 1 TO MM: MX*CK. M)=MZ#CK, M) 

: NEXT M: NEXT K 

3190 ' 

3200 'ROV REDUCTION 

3210 'Input is MX#CK.M) and MPCK) 

3220 'Output is the altered MX^CK.M) and MPCK) 

3230 FOR § = 1 TO MM-1 
3240 FOR K = S+l TO MM 
3250 Q=-MX#CK.S)/MX*CS.S) 

3260 MPCK)=0*MPCS)+MPCK) 

3270 FOR M = 1 TO MM 

3280 MX,*CK.M)=Q*MX#CS.M)+MX*'CK.M) 

3290 NEXT N 
3300 NEXT K 
3310 NEXT S 
3320 ' 

3330 'SOLVE FOR MVCM) 

3340 'Continuing on using MX^CK.M) and MPCK) from before 
3350 FOR S = MM TO 1 STEP -1 
3360 MVCS)=MPCS) 

3370 FOR M = S+l TO MM 

3380 MVCS)=MVCS)-MVCM)*MX*CS. M) 

3390 NEXT M 

3400 M V C S ) = M V C S ) / M X ft (. S . S ) 

3410 NEXT S 
3420 ’ 

3430 ' SOLVE FOR MUCK) 

3440 FOR K = 1 TO MM 

3450 FOR M = 1 TO MM 

3460 M U C K ) =MU (. K ) +M SC K . M ) *M VC M ) 

3470 NEXT M 
3480 NEXT K 

3490 ' PRINT OUT RESULTS 
3500 PRINT "ft" “ V/m V‘ Um " 

3510 FOR M = 1 TO MM 

3520 PRINT M. MVCM). MUCM): PRINT 

3530 NEXT M 

3540 OPEN "0”.#1,AS *This section puts the output onto 

disk. 

3550 PRINT #1,"NN = NN : ” MM = MM 
3555 PRINT tf\ . "RADIUS = '•;A*R; , • radius =”;R 

3560 PRINT P\," " 

3570 PRINT i+l . " . Vm ‘V Um " 

3580 FOR M = 1 TO MM 

3590 PRINT #1. M, MVCM), MUCM): PRINT 

3600 NEXT M 

3610 CLOSE 

3620 END 
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This section contains routines used in the main program. 

1000D ' INVERSION OF A MATRIX 
10010 ' Input is MX£(k. M) 

10020 'Output is MI#CK,M) 

10030 FOR k = 1 TO MM 

10040 MI#(K.K)=1 

10050 NEXT k 

10060 FOR S = 1 TO MM-1 

10070 FOR K = S+l TO MM 

10080 Q=-MX^CK,S)/MX#CS.S) 

10090 FOR M = 1 TO MM ' 

ioioo mx#(k,m)=o*mx*'(s.m)+mx#<:r.m) 

10110 MI^OC. M)=0*MI^(S. M)+MI#(k. M) 

10120 NEXT M 
10130 NEXT k 
10140 NEXT S 

10150 FOR S = MM TO 1 STEP -1 
10160 FOR k = S-l TO 1 STEP -1 
10170 Q=-MX*(k. S)/MX^CS S) 

10180 FOR M=MM TO 1 STE£ -1 
10190 MX^Ck. M)=Q*MX*(S, M)+MX^Ck. M) 

10200 MI*Ck.M)=0*MI*CS.M)+MI,?(k'.M) 

10210 NEXT M 

10220 NEXT k 

10230 NEXT S 

10240 FOR k = 1 TO MM 

10250 FOR M = 1 TO MM 

10260 MI^(k.M)=MI^(k.M)/MX^(k.k) 

10270 NEXT N 
10280 NEXT k 
10290 RETURN 

11000 'MULTIPLY TVO MATRICES 

11010 'Multiply MX-s'Ck.M) times MY^CK, M) in that order! 

11020 'The output is MZ^Ck.M) 

11030 FOR k = 1 TO MM 

11040 FOR M = 1 TO MM11050 MZ^(k.M)=0 

11060 FOR S = 1 TO MM 

11070 MZ*Ck. M)=MZ*(k, M)+MX,VCk. S)*MY#CS Jl) 

11075 IF ABSCMZ^Ck, m 3)<. 0000005 THEN M^Ck.M)=0 

11080 NEXT S 

11090 NEXT M 

11100 NEXT k 

11110 RETURN 

12000 'ADD TVO MATRICIES 
12010 'Add MX^(k.M) to MY^Ck.M) 

12020 'Output is MZ#Ck,M) 

12030 FOR R = 1 TO MM 

12040 FOR M = 1 TO MM 

12050 MZ*Ck.M)=MX*CR,M)+MY#(R,M) 

12060 NEXT N 
12070 NEXT k 
12080 RETURN 
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